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Shearer’s point process and the hard-sphere 
model in one dimension 

Hofer-Temmel Christoph (math@temmel.me)* 


Abstract 

We revisit the smallest non-physical singularity of the hard-sphere 
model in one dimension, also known as Tonks gas. We give an explicit 
expression of the free energy and reduced correlations at negative real 
fugacity and elaborate the nature of the singularity: the free energy is 
right-continuous, but its derivative diverges. We derive these results in 
several novel ways: First, by scaling up the discrete solution. Second, by 
an inductive argument on the partition function a la Dobrushin. Third, 
by a perfect cluster expansion counting the Penrose trees in the Mayer 
expansion perfectly. Fourth, by an explicit construction of Shearer’s point 
process, the unique i?-dependent point process with an i?-hard-core. The 
last connection yields explicit and optimal lower bounds on the avoidance 
function of i?-dependent point processes on the real line. 
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1 Introduction 

The present paper investigates the behaviour of the partition function of the 
one-dimensional hard-sphere model El. also known as Tonks gas [23], at small 
fugacities. Although the one-dimensional case is not as physically relevant as 
the higher-dimensional cases, it is interesting because explicit and beautiful cal¬ 
culations are possible. This paper collects known results and adds new ones to 
complete the picture. 

The partition function has a smallest non-physical singularity at negative 
real fugacity. This fugacity limits the range where high-temperature (equals low 
fugacity) cluster expansions a la Mayer (Ti] are possible and where analyticity 
of the free energy and reduced correlations is easily established. The value of 
the singularity is already known for a long time and has been derived by cluster 
expansion techniques mm- 

We give explicit expressions for the free energy and ratios of the partition 
function at negative real activities between zero and the smallest non-physical 
singularity. Usually, only bounds are given (and possible). Not surprising, the 
Lampert W function mm and two discrete relatives make their appearance. We 
obtain these expressions by different approaches: scaling from the discrete case, 
cluster expansions or alternative viewpoints on the partition function via a com¬ 
binatoric relationship to the avoidance function of a particular one-dependent 
point process. As a consequence, we shed more light on the nature of singular¬ 
ity: the free energy is right-continuous at the singularity, whereas its derivative 
diverges at this location. 

The discrete one-dimensional hard-sphere model lives on Z with a hard-core 
radius parametrised by k. We collect the results about the discrete case EE 
HT] and specialise results from [32]. We want to point out, that in this case a 
purely inductive approach a la Dobrushin [^ yields the full picture and cluster 
expansion is not needed at all. The first way to approach the continuous model 
is to scale the resulting quantities towards a range 1 hard-sphere model on K. 
We recover the ubiquitous Lambert’s W function. The scaling works for the 
critical values, but not for all related quantities. 

The other approaches work directly with the continuous model. In [20] . we 
generalised the inductive approach of Dobrushin [Bj to the continuous setting. A 
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variant thereof, inspired by the calculations in m, identifies the smallest non¬ 
physical singularity directly and yields bounds on some of the related quantities. 

For the continuous case, we extend the tree-operator approach for cluster 
expansions of discrete polymer models from |22j to the hard-sphere model. We 
exhibit a partition scheme of the cluster adapted to the one dimensional setting. 
The singleton trees of the partition scheme have an explicit structure, which al¬ 
lows us to derive cluster expansion expressions exactly, instead of just bounding 
them. This approach was partly motivated by the fact that [8], the current 
best known cluster expansion approach using tree-operators for the hard-sphere 
model, does not attain the critical value in one dimension. We describe the free 
energy at negative real fugacity explicitly up to this singularity. At the fugacity 
approaches the singularity coming from zero, the free energy expression con¬ 
verges, but its derivative diverges. 

Finally, we use the connection between the partition function of a i?-hard- 
sphere model and the avoidance function of Shearer’s point process HBUD]. 
Shearer’s point process is the unique one-independent point process with one- 
hard-core configurations. An alternative proof of a lower bound of the singularity 
comes from an explicit construction of Shearer’s point process. The construction 
is a one-sided variant of a Matern-type deletion rule which seems to seems 
to depend on the chordal (i.e., tree-like) structure of K. See [T^ for the case of 
chordal graphs. As some of the results in the discrete case have been used to 
solve problems in fc-dependent percolation on trees iia, we intend the explicit 
expressions of the continuous case as a tool to treating dependent continuous 
Boolean percolation problems. 

A related area of research are the combinatorics of the virial expansion isms]. 
These are not a topic of discussion here. 

The following subsections introduce the key terms in detail. Section sum¬ 
marises the results. The scaling approach is in section [231 the inductive approach 
in section the tree-operator approach for cluster expansions in section and 
the constructions of Shearer’s point process in theorems and 

1.1 The generating function 

Let (A, d) be a Polish space. A set of points in X is i?-hard-core, if all pairs of 
points have mutual distance at least R. Denote by R-.Q the indicator function 
of i?-hard-core point sets. On a bounded volume A d A and for weight z G C, 
we regard the generating function 

RR{{xi,---,Xn})Wdxi, ( 1 ) 

n=0 

with integration to some standard measure on A. Let A be a Borel set with 
diam(A) < R. There is a deletion-contraction identity 

Z{z, A) = Z{z^ K\ A) + z I Z{z, A\ B{x, R))dx. (2) 

JA 
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For A (B X fixed, the partition function has its smallest root on the negative 
real axis. A key role is played by the quantity 

p* := sup{p > 0 I VA (E A : Z{-p, A) > 0} . (3) 

It is a singularity of the cluster expansion 

OO 

log Z{-p, A) = - ^ 

n—1 

where the Ursell coefficients [53] U{xi,... ,x„) G Nq are combinatorial expres¬ 
sions depending on the whether pairs of points in {xi,...,x„} have distance 
less than R. 

1.2 Shearer’s point process 

A point process ^ on {X,d) is i?-dependent, if, for every pair of subsets Ai, A 2 
with d(Ai, A 2 ) > i?, the projections Ai^ and A 2 ^ are independent. For a given 
intensity, there is at most one i?-dependent point process with an almost-sure 
i?-hard-core. This is Shearer’s point process [iniiinj. 

The avoidance function of Shearer’s point process with intensity p is Z eval¬ 
uated at weight —p. Thus, Shearer’s point process exists, iff it has intensity 
p < p*. li Shearer’s point process exists, it has the minimal avoidance function 
among all i?-dependent point processes of the same intensity. 

1.3 The hard-sphere model 

The hard-sphere model with radius R is the unique Markov point process on 
(A, d) with interaction range R and an almost-sure i?-hard-core. For a bounded 
subset A (E A, fugacity X G [0, oo[ and empty boundary conditions, its partition 
function is Z(A, A). 

The behaviour of log Z(z, A) and ratios of the partition function at different 
volumes in low fugacity (high-temperature) setting are a topic of longstanding 
interest [141117] . Uniform bounds (in both volume and fugacity) in a small disc 
in C around the origin are possible because of the inequality [3 Theorem 2.10]: 

|logZ( 0 , A)| <-logZ(|z|, A). (5) 

Knowledge about the location of p* is crucial to expansions of \ogZ{z^A). Be¬ 
cause p* represents a negative fugacity, it is also called a non-physical singularity. 

1.4 Lambert W function 

The Lambert W function [5| are the solutions of the equation 

W{z) exp(lF(z)) = z 

in the complex plane. For p G [0, j], there is a real branch Wq of solutions in 
[0,1] solving the equation 1F(—p) exp(lU(—p)) = —p uniquely. This solution 


p" f 

— / U{xi, . . .,Xn)dxi 
n\ Ja- 


' dXr^ 


(4) 
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plays a key role in the low fugacity case of the Tonks gas [31 Section 3] and 
appears also in our results. The reason is that the Lambert W function has a 
close relationship with the generating function T of rooted trees on n labelled 
points O Section 2]: T{x) = —W{—x) and T converges, iff x G [0, ^]. 


2 Results 

2.1 The discrete case 

For k G No, we regard the space (Z, |.|) with hard-core radius k + 1 and inte¬ 
gration with respect to the counting measure. The subscript k marks related 
quantities. Instead of point processes, we talk of a Bernoulli random field (short 
BRF), i.e., a collection of {0, l}-valued rvs. This section summaries already 
known results for comparison with the continuous case in section [2.2[ 


We use the convention 0° = 1. Consider the map 


fk : [0,1] 


0, 


X !->■ x(l — x)^ 


For p G 


(/c-k !)'=+! 

^ , let Afc(/o) be the unique pre-image of fk in 


( 6 ) 


’ (/c+l)'=+i 


(k+l) 


k 

Theorem 1. For p G [0, let X := (X„)„gz be a Bernoulli produet 

field with parameter Afc(p). Define the BRF Y := (T„)„gz by 


k 

Yn -.= Xr^Ylil - Xn-^) . 


(7) 


Then Y is Shearer’s BRF on Z, for marginal probability p and distance (fc -k 1). 

7 k 

Shearer’s BRF does not exist for higher marginal probabilities than prpyjfc^. 


The construction is from na, inspired by m and may be seen as a special 
case of the general construction on chordal graphs m- In particular, the BRF 
y is a (fc -k l)-bloek factor of the BRF X. The non-existence follows from (8a) 
in theorem |2j 


We have Zk{p, n) := Z{—p, {!,..., n}), with Zfc(p, 0) = 1 and Zk{p, 1) = 1 —p. 
Let pk(n) be the smallest root of Zk[p, n). If Zk{p, n— I) > 0, then let ak[p, n) := 

. Let Fk{p) := lim ^k{p,n) ^ singularity reduces to 

Zk{p,n 1 ) n—>-oo n 

Pk := sup{p > 0 I Vn G N : afc(p, n) > 0} = inf{pfe(n) | n G No} . 
Theorem 2. For each fc G Nq, we have 


Pk (yt + i)fc+i’ 

Vn G No : Pkin) > Pk{n + 1) with equality if k = 0 , 
lim pk{n) = pI , 
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Vp g] 0, p^.], n S N : ak{p,n) > ak{p,n + 1) with equality if k = 0 , (8d) 

Vp e]0,p;^],n e N : Zfc(p,n) > Zfc(p,«+1), (8e) 

VpG[0,Pfc]: lim afc(p,n) = 1 - Afe(p), (8f) 

n—^-oo 

Vpe[0,p’J]: P’fc(p) = log(l-Afc(p)), (8g) 

lim ^Fk{p) =oo ifkf^O. (8h) 

p-fp*- op 


Statements (8a), (8b) and (^m are already in m- Statement ( |8d[ ) follows 
from the general statement in |21j . with (8e) being a corollary. Statement (|^ is 
from [T3|. Statement (8g) follows from ( |8fp by telescoping and the monotonic¬ 
ity (8d). The proof of statement (8h) is in proposition]^ 


2.2 The continuous case 

The continuous case is on (M, |. |) with hard-core radius 1 and integration with 
respect to the Lebesgue measure. The subscript c marks related quantities. 

Consider the map 


/c : [0, oo[-^ 


0 , 


y ye~ 


(9) 


For p S [0, i], let Ac(p) be the unique pre-image under /c in [0,1]. 

Theorem 3. Let p G [0, -] and ^ be a homogeneous Poisson PP with intensity 
Ac(p). Define a PP p by the density 


p{dx) = 1 £,{dx) = 1 and ^(]a; — 1, a;[) = 0 . 


( 10 ) 


Then p is Shearer’s PP on K with intensity p and distance 1. Shearer’s PP does 
not exist for larger intensities than ^. 

The construction of Shearer’s PP may be seen as a Matern type point pro¬ 
cess [12] with a one-sided simultaneous thinning rule: a point of the Poisson 
point process ^ survives, if it has an empty open unit interval to its left. In 
analogy with the discrete case, we may call p a 1-block factor of the PP f. 


We look at Zc{p,t) ■= Z{—p, [0, t]), with smallest root Pc{r), and ac(p, s, t) := 
z^{p s+t) defined. Let Fcip) '■= hm ^r(p, t) ^ ^ clear from context, 

we omit it. The singularity reduces to 

p* := sup{p > 0 I Vs, t e [0, oo[: adp, s, t) > 0} = inf{pc(t) | t G [0, oo[} . 

Theorem 4. 



Vt G [0, oo[, s > 0 : pc{t) > pc{t s ), 
lim pdt) = p* , 
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Vp £]0, p*], s, s , t, t' G [0, oo[ with s + t' > 0 : 

adp, s, t) > ac(p, s + s' ,t + t'), 

Vp e]0, Pc], t G [0, oo[, s g] 0, oo[: Zc(p, t) > z^p, t + s), 

VpG [0,p*],Vt,sG [0,oo[: ac(p,s,t) > exp(-sAc(p)), 

w ^ rn *1 r logac(p,s,t) , . 

Vp G [0, Pc : hm --- = -Ac(p), 

>-oo 5 

s—vO 


Vp G [0, Pc] : Fc{p) — p 1 + 2pAc(p) + 


P^Ac(p) 


5 

lim —Fc{p) = oo . 
p-s-p*- op 


(lid) 

(lie) 

(l lf) 

(llg) 

(llh) 
(Hi) 


Statements (11b) and (lid) are from m- Statement ( |llh[ ) is from propo¬ 
sition 


sition 


Statement 

Statement 


lli| is from proposition |13| Statement 


llfl is from proposition Ril Statement (11a) follows either 


from propositions and|10| or from proposition and 


llg 


is from propo- 


2.3 Scaling 

This section discusses how the parts of the results for continuous case may be 
seen as scaled versions of the corresponding discrete results. 


The singularities scale: 


{k + l)pfc — ( 1 — 


k i J fe—>oo e 

Even more, the activities and their corresponding inverses, too. If p G [0,1], 
then ^ G [0, ^] and 

k 


{k + l)fk 


y 


k 


= 2 / 1 - 


y 


k -\- \ ) k—yc 


-Gye y = fc{y) ■ 


Hence, for p G [0,1], 




7-^ ^c(p) • 

k—¥oo 


For Shearer’s point process, the scaling is even more obvious: The underlying 
Bernoulli product field scales to the Poisson point process of correct intensity 
Ac(p). The he dependent thinning rule is the same and the length of the free 
interval to the left scales with Therefore, the critical values scale. Note: 
both (1 -I- x)~^ and (1 — cc)^ have the same scaling limit ye~y under x = 


While in the continuous case there is only one fixed point equation, in the 
discrete case there are two: p = x{l — xd is in whereas p = a:(l —= 
()fc(a:) stems from the discrete one-sided tree-operator in section [4^ As both 
(1 -I- and (1 — a;)*’ have the same scaling limit e~y under x = these 

fixed point equations merge in the continuous setup to p = ye' -y (§. 
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2.4 Previous results 

This section lists previous results obtained via cluster expansion. Previous re¬ 
sults in the discrete case are 

• Dobrushin U- pI> ( 2 ifiS+i ~ 2 ^ 

• special-casing Fernandez and Procacci [7]: p* > ~ (71^ 

(maximum of . 

Previous results in the continuous case are 

• Ruelle [T7] p* > l/(2e) 

• special-casing Fernandez, Procacci and Scoppola i: ^ (maxi¬ 
mum of - - - 2 -)- 

l+2/j.+ ^^ 

3 The inductive approach 

3.1 The discrete case 

Let p < • The identity ([^ rewrites as 

o-k(n) = 1-i- - -. 

nLa,(p,n-z) 


We deduce (8d) and show that 


Ofe(p,oo) := lim afc(p,n) 

n—>-oo 


( 12 ) 


exists and fulfils the identity 


afe(p,oo) = 1 - 


(13) 


Ofc(p,oo)'= ■ 

This is just another form of the fixed-point of (§, whence 1 — (p, 00 ) = \k{p)- 

Proposition 5. If k > 1, then 


lim — lim Ofc (p, n) = 00 . 

P^pI~ op n->-cx> 


(14) 


For k = 0, nothing happens. 


Proof. We derive with respect to p in (13), apply it to the result and obtain 

d 1 

— afe(p,oo) = 


j^-ak{p,ooy 


Because of ([^ and ak{Pk, 00 ) = 1 — Xk{ak{p, 00 )) = we have 


lim ^^-«(p,oo)‘ = hC.)*« 


P~^Pk 


V/c + i 


= 0. 




















Hofer-Temmel Christoph 


Shearer’s PP ID 


3.2 The continuous case 

The following proposition is a special case of the continuous Lovasz Local 
Lemma in [20]. The identity ([^ has the form: if s € [0,1[ and t S [0,oo[, 
then 

(15) 

(16) 

(17) 

(18) 


Zc{p, t + s) = Zc{p, t)-p Zcip, (t + x-l) A 0)dx . 

Jo 

If Zc{p,t) > 0, this rewrites into 

adp, s,t) = 1 — p / adp, {1 — x) V tdt + X — 1) A 0)~^dx . 
Jo 


If 


exists, then we have 


adp,s,oo) := lim adp,s,t) 

t —>-00 


adp, s,oo) = 1 — p / ac(p, 1 — a;, oo) ^dx. 

Jo 

Proposition 6. If p < \, then 

Vs, t G [0, oo[: ac{p, s, t) > exp(—sAc(/9)). 

Proof. We omit the trivial case with p = 0. Fix p > 0 and A := Xdp) > 0- Let 
^]^jg jneans that p = Xh'~^ and that diz^ = Xv^dx. We use 


v^-^dx = 


We telescope to reduce to the case s G [0,1[. We proceed by induction over 
[s + tj. We have three base cases and one induction step. 

Case [s + tj = 0: Consider fi{s,t) := 1 — (s + t)p — (1 — tp)v~’^ on [0,oo[^. 
We have /i(0,0) = 0 and 


'Jfi{s,f) = 


-p+ (1 - tp)iy- 


> 


-p + pv 

Therefore, /i is non-negative on [0, oo[^ and 


-7 + (l-7) 

p{v-^ - 1) 


> 


, zds + f) l-{s + t)p 

ac[s^t) = -77^— = -:-^- > Id 




i — tp 


Case [s + ij = l,t < 1: On [0,oop, consider the function / 2 (s,t) := 1 — 
(s -b t)p -b 5 (s -b t - l)2p2 - (1 - tp)i'~‘^ = /i(s,t) + J{.s + t- l)2p^ > /i(s,t). 
Therefore, /2 is non-negative on [0, oo[^ and 

zds + t) l-(s + t)p+J(s + t-ldp^ 
ads,t) = - 77 ^— = -:-^ 




1 — tp 


Case [s + = 1, t > 1: For x G [0, s[, we have t-bx—l<t-bs — 1<1. 

Use (161 and the previous case to get 

ads, t) = 1 — p adl — x,t + X — l)~^dx > 1 — p / d~^dx = v~^ . 

Jo Jo 
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Case [s + tj > 2: Use (161 to get 

ac(s, t) = 1 — p / ac(l — x,t + X — l)~^dx . 
Jo 


For each x G [0, s], we have > [t + s — IJ > [t + a; — IJ. We use (16) 

and the induction hypothesis to obtain 

pl — X 

ac(l - x,t + X - 1) = 1 - p / ac(l - y,t + X + y - 2)~^dy 

Jo 


> 1-p 


pL — X 

/ i'^~^dy = 

Jo 


We plug this into the first expansion and get 

pl — X 


pl — X 

{I — x,t + X — 1) = I — p / 0^(1 — y,t + X + y — 2)~^dy 

Jo 

>1- p [ v^-^dx = . 


□ 


4 Cluster expansion 

4.1 Ursell coefficients 

This section recapitulates how we may express the Ursell coefficients U{xi ,..., a;„) 
as an explicit sum of trees. 

For a finite graph G, let Qo be the spanning subgraphs of G and 7 g be 
the spanning trees of G. A partition scheme S' is a function To ^ Qg, such 
that the intervals {[T, S(T)] :T gTg} partition the poset (Qg,^) into dis¬ 
joint and bounded lattices. The singleton trees of S are its fixed points: Sg '■= 
{TgTg- S{T) = T} 


For Xi,... ,Xn G X, the cluster G{xi, ...,x„) is the graph with vertices [n] 
and edges between vertices i and j with d(xi,Xj) < R. If we have a partition 
scheme of a cluster G{xi, ..., x„), then Penrose’s theorem [TS] implies that 

H spans G(a:i 

Hence, 

OO n Tt 

logZ(-p, A) = - ^ ^ / |SG(xi,....a;„)| ■ 

n=l "'A" 

4.2 The one-sided partition scheme 

This section describes the one-sided partition scheme, a partition scheme adapted 
to a one-dimensional space. Proposition describes its very nice and explicit 
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set of singleton trees. 

Fix xi,... ,Xn & X. Let G := G(xi,..., Xn)- We root G at the vertex 1. Let 
Tn be the set of trees with n vertices. For T G Tn and x G A”, we have 

[T G = [SiT) =T,TC Gixi,Xn)] . 

If X is one-dimensional, i.e., it has a total order compatible with the distance, 
then there is a partition scheme with an explicit description of its singleton 
trees. For T G Tg, let G [n] —No denote the level of the vertices in T, i.e., the 
distance to the root 1. If k > I, let pi be the parent of i. If k > 2, let gi be the 
grandparent of i. 

Proposition 7. There exists a partition scheme One, such that 

n 

[T G Sg] = = 1] + [h > i,Xi< Xp. < XgJ -I- [k > l,Xi> Xp. > 

i=l 

If we regard T in Tn instead oJTg, then we have to add the eluster constraints 
from G: 


[T G 5g] = - Xpi \ < i?] 

i^2 

n 

X = 1] + [h > < Xp, < Xg,] + [k > l,Xi> Xp, > Xg,]^ (19) 

i^2 

We call the partition scheme One the one-sided scheme, because after the 
first-level children have chosen a direction (increasing or decreasing), the subtree 
based on such a child branches only in the given direction, with no additional 
constraint on and between further descendants. 


The one-sided scheme One is an explorative partition scheme |221 Section 
5.2]. An explorative partition scheme has two key ingredients: an exploration 
algorithm (selecting a tree from a cluster) and a tree edge complement partition 
(getting One{T) from a tree T). We state these two ingredients below and omit 
the proof of their correctness and of proposition as the one-sided partition 
scheme is a variation of the returning scheme [22l Section 5.4] and the proofs 
would be slight rewritings of those of [22l Propositions 24 & 25]. The key idea of 
the returning scheme is “The exploration algorithm should select those edges, 
which we want the singleton trees to contain.”. For the one-sided scheme, this 
becomes: Once a direction is chosen, the exploration algorithm below tries to 
go in this direction as far as possible. See the P case in algorithm i If you 
ever have to turn back (i.e., change direction), then the tree edge complement 
partition ensures with (23) that One{T) T. The first level gets a much simpler 
exceptional treatment based on the greedy scheme [5^ Section 5.3]. For the 
remainder of this section, the notation follows [5^ Section 5]. 


Algorithm 8 (Exploration algorithm). Let H G Gg- For every k, let Hk, Tk, 
Uk, Bk and Pk be as in generic exploration algorithm]. The missing parts 
to construct Hk+i from Hk are: 


II 



Hofer-Temmel Christoph 


Shearer’s PP ID 


Call an edge {i,j) G E{C C Pk, Bk) C E{H) progressing (short V), if k = 0, 
or k > 1 and either Xp^ > Xi > xj or Xp- < Xi < Xj. All other edges in 
E{C n Pk,Bk) n E{H) are reversing (short RJ. A vertex is i G Pk progressing, 
if all edges (i,j) are progressing and reversing, if there exists a reversing edge 
{i,j)- Finally we say that a connected component C of Hk\uk progressing, if 
all vertices in C D Pk are progressing and reversing, if C D Pk contains at least 
one reversing vertex. 

For k = 0, we do the same as in the greedy partition scheme \2^ . For k > 1, 
we do as follows: 


If C is an P connected component of Hk\uk- 

(dpb) SELECT C C Sk := C f] Pk- 
(dpi) As C C Ik = ^ REMOVE nothing. 

(dpp) For eachi G CnSk, let ji := argmin{j G Bk : {i,j) G E{Hk)}. SELECT 
ihji)- 

(dpu) For each i G C Cl Sk, REMOVE all (i,j) G E{IIk) with ji ^ j G Bk. 

(dpc) REMOVE all of E{C C Sk) C E{Hk). 

If C is a H connected component of Hk \ Uk ■ 

(drb) SELECT C 0 Sk := {i G C H Pk : i is R}. 

(dri) Cnik = {iGCnPk:i isP}. REMOVE all of E{Bk,iCnIk))nE{Hk). 

(drp) For each i G C D Sk, let ji := argmin{j G Bk : {i,j) G E{Hk) is R}. 
SELECT {i,ji). 

(drur) For eachi G CCiSk, REMOVE every R {i,j) G E{IIk) withji ^ j G Bk. 
(drur) For each i G COSk, REMOVE every P {i,j) G E{CCSk, Bk)r\E{FIk). 

(drc) REMOVE all of E{C C Sk) C E{Hk). 

Algorithm 9 (Tree edge complement partition). Let T G Tg- Let Lk be the k*^ 
level o/T. First we determine if an edge {p{i),i) is a progressing (short P) or 
reversing (short RJ edge: 


d : 


/\{o}\Li^{P,R} 



if C: ^p(i) — ^p(p(i)) — ^p(i) — ^p(p(i)) 

else.. 


( 20 ) 


For k > 2 define the equivalence relation ~(fe) on Lk by 


'i^ik)j ^ d(P{o,i)\{o}) = d{P{o,j)\{o}), (21) 


where the equality on the rhs is taken in {P,R}^“^ between the labels of the 
paths P{o,.) to the root. This implies that an equivalence class consists of either 
only same or only non-same nodes and whence we can extend d to them. For 
completeness let ^(o) be the trivial equivalence relation on Lq and be totally 
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disjoint on Li. The equivalence classes possess a tree structure consistent with 
T; 

* ~(fc+i) j ^ p(^) ~(fc) P(j) 5 (22) 

that is equivalent vertices in L/c+i have equivalent parents in L^. We therefore 
call the parent class o/. 

We partition E\E{T) into AoneC^) WConeClT). Edges incident to the root 
o are treated as in the greedy scheme in i2S^ . 


Let I < k < I, j € Lk, i € Li and e := {i,j) ^ E! \ E{T). Then e G Cone('^)t 
iff one of the mutually exclusive conditions (23) holds: 


[i](/c) ^ ^([o]( 0 ), W(/)), (23a) 

^ > 2 A [j](fc) e D([o](o), [p(p(i))]p_ 2 )) A Xiy^ Xj , (23b) 

^ > 2 A [j](fc) e D([o](o), [p(p(*))]p- 2 )) A Xi= Xj A d{C) = P, (23c) 


where C G P{[j](k), [p(*)](z-i)) the unique class with p(C') = [j](fe), 


^ > 1 A [j](fc) = [p(i)](;_i) A Xj A d{i) = R A j < p(*), (23d) 

^ > 1 A [j](fc) = [p(i)](;_i) A ^ Xj A d{i) = P, (23e) 

^ > 1 A [j](fc) = [p(i)](;_i) A x^= Xj A d{i) = P A j < p(i). (23f) 

And e G AoneCT^), iff one of the mutually exclusive conditions ( |24[ ) holds: 

[j]{k) = \i]{i) , (24a) 

l>2 A [j](fe) G P([o](o), [p(p(*))](i- 2 )) A Xi= Xj A d{C) = R, (24b) 


where C G P([j](fc)) [p(*)](;-!)) the unique class with p(C') = [j](k)> 


1>1 A [j](fe) = p([i](i)) A Xif^ Xj A d{i) = R A j > p{i), (24c) 

^ > 1 A [j](fe) = p([*](i)) A Xi = Xj A d{i) = R, (24d) 

^ > 1 A [j](fe) = p([*](i)) A Xi = Xj A d{i) = P A j > p(i). (24e) 


4.3 Tree-operator series 

This section builds heavily on [52 Section 4], in particular [52 Proposition 7]. 
We encode the properties of the singleton trees of the one-sided partition scheme 
in proposition [^ into convergence conditions for cluster expansion series. In this 
section, we often give three expressions: a generic one without subscript, fol¬ 
lowed by the specialised ones for the discrete and continuous cases. 


In the tree-operator of the Penrose partition scheme Pen [15] [22 Sections 
5.1 & 5.3] a vertex has s children with weight G{s): 


GIs) := 


[s = n] 


n=0 J]-R,Rl’- 

Gk{s) '■= [s = 0] -|- {2k -|- 1)[s = 1] -p 


Py{d)y, 


k + 1 
2 


= 2] 
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Gc{s) — [s — 0] + 2[s — 1] + 
leading to a tree-operator defined by 

OO 


= 2 ] 


s=0 


gki/J-) 1 -b (2fc + l)fj, + 


9c{g) 1 -b 2/i -b 




fc -b 1 
2 

2 


The function 


0(Ai) := ^ Qk{g) ■= —^ 
9KN 9k W 

has a global maximum in [0,oo[ at 

2 ^ \/k{k -b 1) 


flc(M) := 




9c{9) 




k{k + l)j ^2k{k + l) + 2k + l 2-bV2’ 


In the tree-operator of the one-sided partition scheme One a non-root vertex 
has s children with weight H (s): 


OO r 1 

n\ 


n=0 


'[O.-Rp 


dV: 


E 1 = E 1 = 


n—Q 


*=1 a;ie[fc-|-l]\{a:i,...,a;i_i} + 


k -\-l 


[s = nj 
d 


n! 




Hc(s) := 

n—0 

leading to a tree-operator defined by 

OO 

d(fJ-) := 


dy = 


Si 


5 = 0 

fc+1 


OO K-Pi /7 I -j \ 

hUM):=E^^(^)d^ = E{ )m^ = (1+m)^+' 

5=0 5=0 ^ ^ 

^ ^ TJS..S 

Mm) :=E^^W^'=E^ 


5=0 


5=0 


The function 


t)(M) := 


M 


f)fc(M) := 


hkiy) 


= e^i^. 


f)c(M) := 


M 


hciy) 
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has a global maximum in [0,oo[ at 

Consider the fix-point equation 

^j. = ph{n) p. = phk{p) p = phc{p). 

As the derivative of the rhs is positive for p < ^ and p < 1: 

dp p dp 

dp 1 — kp dp ^ — P 

Therefore, the inverse p{p) is well-defined up to the maximum. In the continuous 
case, the fix-point equation is p = pe~^ = fdp), the derivative is|^ = {l—p)e~^ 
and the inverse is pdp) = ^dp)- 

Consider the operator Tp ■. p^ ph{p). By [22l Proposition 7], if p < p*, then 
the generating function Q of rooted, one-sided trees converges and equals 

OQ Tl 

Q(p) := lim r"(0) = XI ^ XI 11 ' 

n—>-oo ^ Tl' 

n=l TeTn i^l 

For n S No, let 

n 

D{p,n):= X G(si)n^r(s.)- 
Ter„ i=2 

The full series P of all singleton trees for the directed scheme needs and extra 
step for the first level (children of the root): 

00^1 n 

n—l T^T'n i =2 n —1 

We also have the truncated series 

N n 

Pn{p) ■■= X ■ 

n\ 

n—l 

Due to measurability, the resulting structure differs slightly between the 
discrete and continuous cases. In the discrete case, we have 

Qk{p) = ppk{p) 


and the full series Pk fulfils 

-Pfe(p) = P + Dfe(p)-I-2fcQfe(p)-I- ^ • 

The rhs is not gk{Qk(,p)), because a child of the root staying at the same integer 
than its parent has not chosen a side yet. Thus, 

Pk{p) = p ‘^^Qk{p) ^ 2 ) Qk{p)^'^ ■ 
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In the continuous case, we have 

Qc{p) = PPc{p) = pK{p) 

and the full series Pc is 

Pc{p) = pgdQcip)) = P + 2Qc(p) + . 

Proposition 10. If p > then Qdp) and Pc{p) diverge. 

Proof. The trajectory of 0 under the map ^ i—> pe^ diverges. □ 

4.4 The shape of the free energy 

Proposition 11. For p € [0, ^], we have 

Fc{p) := lim _ ^ (25) 

t—fCXD t 

Proof. Fix t > 0 and p S [0, j]. For T G 7)i, let 

Td) ■= [ [T ^ Sa(s)]dx. 

We also have 

R{t,n) := [ |5G(s)|df= [ X! & SG{x)]dx = ^ T(t) 

and 

OO ^ n OO Yl 

S{p,t) := - log Zcip,t) = / \SG(s)\dx = '^^R{t,n). 

Summing in from the leaves a la Cammarotta |3] we bound 

n 

T{t) < tGc{sl)Y[Hc{s^) and R{t,n) < tDc{p,n). 

i^2 

This implies 

Ti 

S{p,t) <tY] ^DciP, n) = tPcip). 

n\ 

n—l 

On the other hand, if the depth of the tree T is at most d and the root label 
xi G [d — 1, r — d + 1], then the vertex labels of T are contained in [0, t]. Hence, 

n 

T{t) >{t- 2{d - 1))+Ge(si) J] Hds,). 

i=2 

A size n tree has at most depth (n — 1): 

Rit, n) > (t — 2{d — l))~'~[n < d]Dc{p^ n). 
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Hence, 


n 

S{p,r) > V —{t-2{d- l)) + [n < d]Dc{p,n) 
n! 

n—\ 

^ Tl 

= it- 2{d - 1))+ ^ ^Dc{p, n) 

n—1 

= (t-2(d-l))+Pd.c(p). 

For small £ > 0, we may choose first d^ such that Pd^^dp) > (1 — d^dp) and 
then tg such that (tg — 2{de — ^ (1 ~ ^)^e- Thus, for all d> d^ and t > tg, 

we have 

Sip, t)>it- 2{d - 1))+Pd,c(p) > (1 - sftPdp) ■ 

□ 


Proposition 12. For p G [0, ^], we have 


lim - 

s—>-oo s 

Proof. Fix p g] 0, i], as the statement is trivial for p = 0. We have 
- log ads, t) = - log zds + t) + log zdt) 


= y- 

J[o,s+t]\[o,s] 

OO 

= 1:'’ 


n\ 


= P 


0 v/[l/,s + t]n--l 

P" f 


/ E 

lo ^0 ''■ 


|5G(a-)|dx 

\SG{y,x)\dxdy 
\SGiy,s)\dxdy 


By Proposition 7], we have that 


n—0 


n+1 


Ac(p) 


E „l / \^Giy,x)\dx-'^ ^ i7(sj) _ 

'[0,oo[”’ ^—n ’ — 1 > 


n\ 

TL — 0 T^'T'n-\-l i — ^ 


On the one hand, we majorate as follows: 


-\ogads,t)<p 


'0 n=0 

P_ 
n\ 


'[y,° 


\SG{y,x)\dxdy 


^ p ' V j 

= ^®E;t/ \SG{y,x)\dx 


n—0 


'[ 0 ,' 


= sXdp) ■ 

On the other hand, we minorate as follows: 

PJ 


\ ' P" f 

-\ogads,t)>p E “T / \SG{y,x)\dxdy 

do d[y,oo[- 
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pU r 

= \SGiy,x)\dxdy. 

„=0 


Therefore, the limits 


lim lim = ,i„, MpI = xAp) 

s—^0 1—>-oo 5 s—^0 S 


and 


logac(s,t) p" 

lim lim-= lim p > —- 

s t^o n\ 

n—Q 


t—^oo s—>-0 


\SG{y,x)\dxdy = \c{p) 


are equal. 

Proposition 13. ITe have 

d 

lim ^Fc{p) = oo . 
P-^P*- op 


Proof. We have 

This explodes at as Ac(^) = 1. 


Q>'a{p) 

1 - Ac(p) ■ 


□ 


□ 
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